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The termination structures of the jets of Fanaroff & Riley (FR) galaxies are observed 
to produce extended non-thermal emission in a wide frequency range. The study of these 
structures can provide valuable insights on the conditions for particle acceleration and 
radiation at the shock fronts. We have studied the thermal and non-thermal emission that 
can be expected from the jet termination regions of Fanaroff & Riley type I sources. The 
broadband emssion from these galaxies has been recently extended to include the high- 
energy gamma-ray domain, owing to the Fermi detection of Cen A lobes. Exploring the 
physics behind the jet/medium interactions in FRI can provide valuable insights on the 
conditions for particle acceleration and radiation in the jet termination shocks. Making 
use of the results of a fully relativistic numerical simulation code of the evolution of a 
FRI jet we model the expected radiative output and predict spectra and lightcurves of 
both thermal and non-thermal emission at different source ages. 
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1. Introduction 

Fanaroff-Riley sources of type I and II (FRI-II) are radio loud Active Galactic Nu- 
clei (AGN) displaying large-scale jets that interact with the surrounding ISM/ICM 
mediumffl The study of the extended emission in FRI/II sources can be used to 
derive the properties of the flows and the surrounding medium, and to determine 
the conditions for particle acceleration. In this regard a simple evolutionary model 
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can be used to conclude that thermal X-ray emission produced in the shell could be 
detected only for dense enough media even for powerful FRII jetspl whilst thermal 
MeV radiation from the cocoon region would be significant only for sources with 
ages <C 10 7 yrsffl Non-thermal radiation can also be produced in the jet/medium 
interaction regions. X-ray synchrotron and/or IC emission has been observed in 
several sources (see e.g. |H [5]), and the Fermi Collaboration recently reported the 
detection of extended GeV emission in the radio lobes of Cen A, demonstrating 
that acceleration up to very high energies is indeed taking place in the disrupted 
jet region.^ Numerical simulations, on the other hand, have demonstrated to be a 
powerful tool to study the jet/medium interactions. Perucho & Mar ti (2007) per- 
formed a simulation aimed to test the FRI jet evolution paradigmPEl We computed 
the thermal and non-thermal emission from FRI sources obtained with a radiative 
model adapted from Bordas et al. (2009) that takes as input the results of those 
simulations as input and the radiation model used to characterize the conditions 
of the shocked regions!^ Below, we briefly introduce the numerical simulations and 
the radiation model used. We then discuss the obtained results and the relevance 
of the thermal and non-thermal radiation at different energy bands. For a more 
extended discussion, the reader is referred to Perucho & Marti (2007) and Bordas 
et al. (2011)!3 ? 

2. Hydrodynamical simulations and radiation model 

The results of the numerical simulations of Perucho & Marti (2007) have been 
used to characterize the physical properties of the jet/medium shocked regions.^ 
The jet is injected in the numerical grid at 500 pc from the active nucleus, with a 
radius of 60 pc. The ambient medium, composed by neutral hydrogen, has a profile 
in pressure (see, e.g.. I13|) . The jet, leptonic in composition, has an initial velocity 
Vjq = 0.87c, a density ratio with respect to the ambient Pjo/p&o = 10 -5 , a pressure 
ratio with the ambient Pjo/PaO — 8, and temperature 4 x 10 9 K, resulting in a 
kinetic luminosity Lj = 10 44 erg s _1 . Further details on the simulation code can be 
found in [71 Regarding the radiative code, we used a simplified one-zone model to 
study the emission properties of both shell and cocoon regions. The injected non- 
thermal luminosity is 10% of the total jet kinetic luminosity Lj. The magnetic field 
B has been fixed taking the magnetic energy density to be 10% of the ram/thermal 
pressure. Concerning particle acceleration, the recollimation shock is assumed to be 
the accelerator of particles in the cocoon region. We use a relativistic approach in 
this case, with an acceleration rate E = r/qBc with rj = 0.1. For the bow shock, we 
have adopted instead a non-relativistic approach in view of the lower shock velocities 
~ (1 — 2) x 10 8 cm s _1 . The particle energy distribution at a given time, N(E, t STC ), 
is calculated considering the injection (oc E~ p , with p — 2.1), at any time t STC . 
Maximum particle energies, E max (t), vary also along time, through the magnetic 
field, the accelerator size and the shock velocity dependence. Further details on the 
radiative code can be found in [9] and [TUl 
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3. Results 

3.1. Non-thermal emission 

The non-thermal spectral energy distributions (SED) for the cocoon and the shell, 
at t src = 10 5 , 3 x 10 6 and 10 8 yr, are shown in the left and center panels of Fig.Q] 
The obtained radio and X-ray synchrotron luminosities in both regions are at the 
level of 2 x 10 41 erg s _1 . The approximate constancy of the X-ray luminosities with 
time is due to the fact that particles have reached the steady state at t STC through 
synchrotron cooling and the assumed constancy of L nt . The decrease of B with time, 
and therefore the growth of t S yn, is compensated by the increase of time available 
for cooling. The IC luminosity grows as long as this process becomes more efficient 
compared to synchrotron and adiabatic cooling, which is shown by the decrease 
of us/ttrad from « 5 x 10 3 (10 5 yr) to 4 (10 s yr). The cocoon and the shell have 
similar HE luminosities, but the cocoon is a few times brighter at VHE than the 
shell due to the higher maximum energies that particles can attain in the former. 
In both regions the bolometric IC luminosities grow similarly with time, reaching 
~ 10 42 and 10 41 erg s _1 at HE and VHE, respectively. The lightcurves for the 
luminosities in radio (at 5 GHz), X-rays (1-10 keV), HE (0.1-100 GeV) and VHE 
(0.1-100 TeV), for the cocoon and the shell, are presented in Fig. [2j They show in 
more detail the time behavior of the non-thermal radiation at different wavelengths 
discussed above. The complex and smooth shape of the lightcurves, most clear for 
the HE and the VHE emission, is a consequence of the hydrodynamical evolution of 
the whole interaction structure propagating in an inhomogeneous external medium. 



3.2. Thermal emission 

Thermal emission has been studied for the shell, much denser than the cocoon. 
We have only considered thermal Bremsstrahlung emission, although we note 
that recombination radiation could also be significant (with luminosities L ICC > 
10~ 21 n c n mn T~ x l 2 V s h Z 4 erg s -1 , where t^h is the shell volume filled with elec- 
trons and ions with particle densities n c and rij on , respectively, and eZ is the ion 
charge^!). We considered two separate regions: one, cooler (ultraviolet -UV-/soft X- 
rays) but brighter, corresponds to the averaged shell conditions (Av), and another 
one, fainter but hotter (hard X-rays), corresponds to a region close to the apex 
of the bow shock, with properties close to those given by the Rankine-Hugoniot 
(RH) conditions. The volume of the latter region is about 3 — 4% of that of the 
whole shell. Figure [TJ right panel, shows the thermal SEDs for t SIC = 10 5 , 3 x 10 6 
and 10 s yr. The slowdown of the bow shock and the velocity dependence of the 
postshock temperature, cx u 2 s , leads to a decrease in the peak of the thermal emis- 
sion with time. However, the overall luminosity increases as the source gets older, 
from 10 39 erg s _1 to few times 10 erg s _1 , since the Bremsstrahlung time-scale 
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Fig. 1. SEDs of the synchrotron and IC (left) and Bremsstrahlung emission (right) for 
t stc = 10 5 (solid line), 3 x 10 6 (long-dashed line) and 10 s yr (dotted line). Both shell and 
cocoon contribute to the overall non-thermal emission. Thermal radiation is computed only 
in the shell, and has been splitted into one component corresponding to the average shell 
properties (Av, black thick), and another one produced in the shell region with conditions 
similar to those of Rankine-Hugoniot (RH , green thin). Only one jet/medium interaction 
region is shown. 



is iBrcmss = 10 3 T 1 / 2 n~ 1 >> t SIC at any source age, with the shell and the hot 
postshock region components peaking from soft X-rays to UV and from hard to soft 
X-rays, respectively. Thermal Bremsstrahlung increases from t src = 10 5 to ~ 10 6 yr, 
then the luminosity slightly decreases until i src ~ 3 x 10 6 yr, when a transition in the 
external medium from the denser galaxy core to the rarefied galaxy group medium 
is produced, and increases again afterwards. The component Av dominates the 
thermal bolometric luminosity in young sources, but the component RH becomes 
similarly bright at i src ~ 10 8 yr. 



4. Discussion 

Our results can be discussed on the ground of the thermal and non-thermal fluxes 
predicted at different energy bands. In radio, the cocoon dominates with fluxes up 
to ~ few x 10" 12 (d/100 Mpc) -2 erg cm -2 s" 1 ~ 10 Jy at 5 GHz from a region 
of few times 10'(<i/100 Mpc)" 1 of angular siz e. Th is radio fluxes are similar to 
those observed for instance in 3C 31 and 3C 15 The shell radio flux is sligthly 
below the cocoon one, although limb brightening may enhance its detectability. The 
radio lightcurve is quite steady, since the accumulation of electrons compensates the 
weakening of the magnetic field with t src . 

At X-rays, the different thermal emitting regions in the shell would lead 
to a bolometric flux from ~ 10" 15 (d/100 Mpc)" 2 (10 5 yr) to a few times 
10" 13 (d/lQO Mpc)" 2 erg cm" 2 s" 1 (10 8 yr) in 3C 15. The hard X-rays would 
be located at the apex of the bow shock, with a typical angular size of 
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Fig. 2. Computed non-thermal lightcurves of the radio (5 GHz x L5 GHz> top-left), X- 
ray (bolometric: 1-10 keV, top-right) and gamma-ray emission (bolometric: 0.1-100 GeV, 
bottom-left; bolometric: > 100 GeV: bottom-right) in the age range t src = 10 5 — 10 8 yr. 

a few V (d/100 Mpc) -1 (for t SIC ~ 10 8 yr), whereas lower energy X-rays 
would come from the whole shell, with of a few 10' (d/100 Mpc) -1 . Regard- 
ing non-thermal X-rays, the cocoon dominates the total output, with fluxes ~ 
10~ 13 (d/100 Mpc)" 2 erg cm -2 s -1 , although limb brightening effects may again 
increase the shell brightness. In fact, in the case of Cen A, the shell seems to 
be the dominant source of non-thermal X-raysP This difference could be ex- 
plained by a higher E max in the shell of that source. The non-thermal X-ray fluxes 
~ 10~ 14 erg cm -2 s -1 of 3C 15 (see e.g. |2]) imply a non-thermal luminosity of 
« 1.3 x 10 41 erg s -1 at 300 Mpc, in good agreement with the values predicted here. 
Furthermore, synchrotron emission concentrated around the recollimation shock is 
compatible with the large-scale jet X-ray emission found in 3C SIM^ 1 If a strong 
recollimation shock is indeed the origin of these large-scale jet X-rays, then the 
hypothesis that jet disruption in 3C 31 is caused by shock triggered instabilities 
is favored against stellar wind mass-load!^ Finally, we note that for sources older 
than those considered here and/or lower £?-values, the synchrotron emission would 
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be less relevant and IC would dominate the X-ray output. 

The moderate velocities of the bow shock found in the numerical simulations 
prevents thermal photons from reaching energies as high as those discussed in [3] 
In the gamma-ray domain, the SED is close to flat at HE, and becomes steeper 
at VHE. We have not accounted for the gamma-ray absorption induced by the 
extragalactic background light, which would become significant at distances larger 
than 100 Mpc. The overall emission increases with time mainly due to the increasing 
efficiency of the IC channel (considering the cosmic microwave background only) as 
compared with synchrotron and adiabatic losses. Gamma-ray fluxes for a source 
with i src ~ 10 8 yr are around ~ 10~ 12 (d/100 Mpc) erg cm" 2 s _1 . At HE, such 
a source may require very long exposures to be detected by, e.g., Fermi, unless it 
is very nearby (as it is the case of Cen A, see |6]) or shows a higher non-thermal 
efficiencies and/or jet power. At VHE, the fluxes could be observed by the current 
Cherenkov telescopes, although the extension of the source, of tens of arcminutes at 
100 Mpc, and the steep spectrum above ~ 100 GeV, may make a detection possible 
only after a long exposure. The forthcoming CTA observatory may however allow 
the detection of VHE emission from FRI jet lobes, and possibly carry out energy- 
dependent morphological studies. 

From our study we conclude that, for moderate non-thermal luminosities, radio 
lobes of FRI radio galaxies are good candidates to be detected in the whole spectral 
range, with the radiation appearing extended in most of the energy bands. Soft 
X-rays will be likely dominated by synchrotron emission up to ages ~ 10 s yr, with 
IC tending to be dominant for older sources. Thermal X-rays seem unavoidable 
and may dominate in hard X-rays in old sources even if a non-thermal component 
is present. The low surface brightness may require long observation times for the 
detection in X- and gamma-rays, although the steady nature of these sources can 
help in this regard. 
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